Infrasonic Detection of a Large Bolide over South Sulawesi, Indonesia on October 8, 
2009: Preliminary Results 

E. A. Silber • A. Le Pichon • P. G. Brown 


Abstract In the morning hours of October 8, 2009, a bright object entered Earth's atmosphere over 
South Sulawesi, Indonesia. This bolide disintegrated above the ground, generating stratospheric 
infrasound returns that were detected by infrasonic stations of the global International Monitoring 
System (IMS) Network of the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO) at 
distances up to 17 500 km. Here we present instrumental recordings and preliminary results of this 
extraordinary event. Using the infrasonic period-yield relations, originally derived for atmospheric 
nuclear detonations, we find the most probable source energy for this bolide to be 70 ± 20 kt TNT 
equivalent explosive yield. A unique aspect of this event is the fact that it was apparently detected by 
infrasound only. Global events of such magnitude are expected only once per decade and can be utilized 
to calibrate infrasonic location and propagation tools on a global scale, and to evaluate energy yield 
formula, and event timing. 
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1 Introduction 

Medium sized Near Earth Objects (NEOs) (>10 m diameter) may penetrate deep into the atmosphere, 
though rarely, and cause significant damage on the ground (Chapman, 2008) and could potentially 
perturb regional climate trends (Toon et al. 1997). However, currently available models cannot 
accurately define the critical impactor size at which the regional climate is affected (Bland and 
Artemieva, 2003). A part of the problem is limited observational data, as records of significant NEOs 
are scarce. Therefore, various observational methods, including infrasound, are critical to re-evaluate 
airburst models and determine with more accuracy the size at which an object can influence the local 
climate. 

Records of significant NEO impacts are rare. Klekociuk et al. (2005) and Arrowsmith et al. 
(2008) report multi-instrumental observations of two different impactors with energies of 20-30 kilotons 
of TNT (1 kT = 4.185 x 10 12 J) occurring in the fall of 2004, while Brown et al. (2002) present infrasound 
data for two somewhat less energetic events over the Pacific in 2001. In all cases these events occurred 
over open ocean and much of the energetics information was compiled from records of the associated 
airwaves detected by infrasonic stations. 
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Infrasound is low frequency sound extending below the 20 Hz hearing threshold of the human 
ear and just above the natural oscillation frequency of the atmosphere (>0.01 Hz, Brunt-Vaisala 
frequency). It has the ability to propagate over long distances with very little attenuation, thus enabling 
the study of remote explosive sources (Hedlin et ah, 2002). The International Monitoring System (IMS), 
operated by the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO), features as one of its 
monitoring technologies, a global network of 42 fully certified infrasonic stations designed to detect 
nuclear explosions (CTBTO web: http://www.ctbto.org). 

Bright meteors (also known as fireballs) fall into the category of events that can be detected and 
consequently studied using infrasound (ReVelle, 1976, 1997; Brown et al., 2002a). Fireballs are 
produced by large meteoroids which may penetrate deep into the atmosphere and generate a cylindrical 
blast wave during their hypersonic passage, which decays to low frequency infrasonic waves that 
propagate over great distances (ReVelle, 1976; Edwards, 2010; Le Pichon et al., 2002a, Brown et al., 
2002; Brown et al., 2003). Global impacts detected infrasonically can provide a valuable tool in the 
estimation and validation of the influx rate of meter sized and larger meteoroids (Brown et al. 2002; 
Silber et al. 2009). Very often, infrasound offers the only available record when it comes to major 
impacts over open ocean. Infrasound observations can provide crucial trajectory and energetics 
information for interesting events which otherwise lack such information (e.g. the Carancas crater 
forming impact in Peru in 2007 (Brown et al., 2008; Le Pichon et al., 2008)). Here we present evidence 
that a significant NEO impact occurred on October 8, 2009 over South Sulawesi, Indonesia based 
primarily on infrasonic recordings of the blast wave detected across the globe; this may have been one 
of the most energetic impactors to collide with the Earth in recent history. 


2 The Indonesian Bolide 

At 2:57 UTC (10:57 a.m. local time) on October 8, 2009 a loud rumbling sound and ground shaking 
startled the people of the town of Bone, South Sulawesi, Indonesia (4.5°S, 120°E). Eyewitnesses who 
ran out their homes in fright saw a very bright object flying across the sky, subsequently disintegrating 
in the mid air, leaving a thick dusty smoke trail behind (Surya news report, in Indonesian: 
http://www.surya.co.id/2009/10/09/ledakan-misterius-guncang-sulsel.html). A news article stated that 
there are reports from local residents that the surviving remnants of the object may have crashed into the 
sea (Surya news report, in Indonesian: http://www.surya.co.id/2009/10/09/leda kan-misterius-guncang- 
sulsel.html). 

Shortly thereafter, the national media, including Metro TV of Jakarta and two news agencies, 
The Jakarta Globe and The Jakarta Post, released a number of reports, including an amateur video of the 
smoke trail (Figure 1). Features and the appearance of the smoke trail are consistent with dust trails of 
other fireballs observed in a similar manner (e.g. the Tagish Lake fireball (Hilderbrand et al, 2006)), 
indicating a probable meteoritic origin of the event. As per The Jakarta Globe, the airburst caused 
damage to several houses in Panyula village (The Jakarta Globe, available at: http:// 

www.thejakartaglobe.com/home/astronomer-sulawesi-blast-bigger-than-atom-bomb-and-caused-by-met 
eorite/338073) and the police department in Bone was flooded with reports of audible sounds extending 
as far as 11 km from Latteko, Bone district, South Sulawesi (The Jakarta Globe: available at: 
http://thejakartaglobe.com/home/mysterious-explosion-panics-locals-in-south-sulawesi-police-still- 
investigating/334246). Unfortunately, there was one casualty, a 9 year old girl with an underlying heart 
condition who went into cardiac arrest upon hearing the thunderous sounds (The Jakarta Globe, 
available at: http://www.thejakartaglobe.com/home/astronomer-sulawesi-blast-bigger-than-atom-bomb- 
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and-caused-by-meteorite/338073). Initially, local people speculated that the event was caused by a 
falling airplane; however, South Sulawesi Police spokesman Sr. Comr. Hery Subiansauri confirmed that 
no aircraft was involved nor any other air incident had occurred. The extraterrestrial nature of the event 
was confirmed by Thomas Djamaluddin, head of the Lapan Center for Climate and Atmosphere Science 
(The Jakarta Post, available at: http://www.thejakartapost.eom/news/2009/10/08/blast-may-be-result- 
falling-space-waste-or-meteorite-lapan.html). 

Upon scrutinizing scrutinizing these reports, we undertook a thorough investigation of infrasonic 
records of all IMS infrasound stations to search for possible signals from the air explosion. 


Figure 1 . A screenshot from Metro TV news report showing an amateur video of the smoke trail, twisted by the 
wind (You Tube, available at: http://www. youtube. com/watch?v=yeQBzTkJNhs&videos=jkRJgbXY-90). 


3 Data Processing and Analysis 

We were able to examine a total 3 1 infrasound stations in the IMS network which were providing data at 
the time of the event. Probable signals originating from 4.5°S, 120°E were detected at 17 IMS stations 
(Figure 2), which we correlated with the event. Table 1 summarizes data from all stations which 
detected the signal, sorted by distance. The signal was extraordinary in two aspects: first, it was detected 
by many infrasound stations, some of which are at extreme distances (>17,000 km), and second, that 
most of the signal energy is contained in very low frequencies, indicative of a source yielding very high 
energy. Infrasonic signals were analyzed using two independent methods, Matseis 1.7 (Harris and 
Young, 1997; Young et ah, 2002) and Progressive Multi-Channel Correlation Method (PMCC) (Cansi, 
1995). 

First, infrasound data across each station have been array processed in windows (typically of 30- 
60 second length) to search for coherent signals with consistent back-azimuth measurements for several 
adjacent windows using the analysis package Matseis 1.7 (Harris and Young, 1997; Young et ah, 2002). 
To determine the arrival azimuth for a coherent signal, we used the standard method of cross-correlating 
the output between each sensor of an array and performed beamforming of the signals across the array 
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Figure 2. A global map (courtesy of CTBTO, web: http://www.ctbto.org) showing all stations (black circles) that 
detected the Indonesian bolide event circled in red. 


Table 1 . Summary of all detections, sorted by distance. We include results for two methods of signal detection 
(MatSeis and PMCC). 
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I39PW 
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(Evers and Haak, 2001). A sample output is shown in Figure 3. In total 15 positive detections were 
identified in this way, using the approximate location and timing from media reports and expected 
typical stratospheric propagation speeds as a guide to isolate the period of most probable signal arrival 
on each array. This procedure was repeated for multiple bandpasses to try and isolate any coherent 
signal from the station noise. 


so 

60 

40 

20 

0 

1 

0.5 

0 

360 

270 

ISO 

90 

0 

1 

0.5 

0 

- 0.5 

-1 


07 : 50:00 08 : 00:00 08 : 10:00 08 : 20:00 08 : 30:00 08 : 40:00 

2009 / 10/09 


Window Parameters 

Pass Band Parameters 

Slowness Parameters 

Windowed Values & S.D. 

Duration (s) 

30.0 

Low frequency (Hz 

0.05 

Max. Slown. 

400.0 

Corr. 

0.887 i 0.059 

Overlap (%) 

50.0 

High frequency 

1 

#ofSlown. 

40 

Fstat 

40.546 30.335 

Number of 

240 

Order (integer) 

2 



Vel. 

0.336 !; 0.000 


196.189; 0.000 

| Animate 1 I SendtoMap | 

[Calculate 1 0 | Write File | C:\Program Files1MATLAB71\work1ai! I Cancel | 


V .’•>* * . \*%,*V* * *' •• 



♦ **♦.**♦.♦ **** * / *.♦ *•* * * / ♦ ** *** ****** ******* 
1 1 1 

i i i 

»*s***++*#S ******* ♦* *♦**♦* ^ + ♦ ****** ******** **** * *V* 

1 1 1 1 


, 1 1 1 

* * ** t \ * 

♦ ♦♦ ,** * * * * 

1 ! 1 l *♦ * l 

1 1 * * 1 - 
m'tWHV.Vv.’V,,' ******* ** ***** «♦ * ♦ 

1 * 1 


MatSeis Upds 
1 1 

dej 

i i 

I45H1BDF 

1 1 1 


Figure 3. An example of signal observed at I45RU, located 5543 km from the source. The top window is the F- 
statistic, a measure of the relative coherency of the signal across the array elements in any particular window, the 
second window represents the apparent trace velocity of the acoustic signal across the array in the direction of the 
peak F-stat, while the third window shows the best estimate for the signal back-azimuth in the direction of 
maximum F-stat for each window. The fourth window shows the bandpassed raw pressure signal for one array 
element. 
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The second method, PMCC, for analysing the data, sensitive to coherent signals with very low 
signal-to-noise ratio (SNR), yielded positive detections at a total of 16 IMS stations. This technique has 
been successfully implemented in detections of other bolides (cf. Arrowsmith et ah, 2008), as it searches 
for coherent signals in both frequency and time windows, selecting detections of similar parameters to 
identify coherent signals (Brachet et ah, 2010) (Figure 4). 



Figure 4. Results from array processing using the PMCC algorithm for the IMS station 145RU. The top window 
gives the observed azimuth, while the middle window represents the trace velocity of the signal. The bottom 
window shows the bandpassed raw pressure signal for one array element. 


We have also established a geolocation using the nine closest stations (Figure 5) by utilizing a 
non-linear system of equations describing the propagation of the detection waves through the 
atmosphere, where the inverse location algorithm is based on Geiger's approach (1910). The location 
results are obtained assuming a homogeneous half-space with a typical celerity value of 290 m/s for 
each individual phase without azimuthal correction (Brown et ah, 2002). In order to detennine the 
location errors, the 95% confidence ellipses are estimated by repeatedly running the linearized least- 
squares inversion with arbitrary sub-sets of the input data within ±10° and ±30 m/s ranges of 
uncertainties for the azimuths and celerity, respectively. 

The maximum peak-to-peak amplitude was detennined by bandpassing the stacked, raw 
waveform using a second-order Butterworth filter and then applying the Hilbert Transform (Dziewonski 
and Hales, 1972) to obtain the peak of the envelope. The filter cutoff frequencies were typically 0.05 Hz 
for the low frequency and up to 2.1 Hz for the high frequency (with few exceptions) and were 
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determined using a power spectral density (PSD) method where the signal segment of the waveform was 
superimposed over the average of the prior and post background noise (of equal length), all being 
divided into equal windows (50-170 seconds in length, depending on station), establishing a frequency 
band which lies above the noise. Therefore, the low and high frequency cutoffs would be selected where 
the signal rises above the noise on the low end or descends into the noise on the high end of the 
spectrum, respectively. 



Figure 5. Map showing the geolocation. The best fit solution was obtained using nine stations closest to the 
Indonesian bolide event. 


To measure the dominant period at maximum peak-to-peak amplitude, two independent 
techniques were employed. First, the dominant period at maximum frequency was acquired from the 
residual power spectral density (PSD) obtained using the method described above, except the noise PSD 
was subtracted from the signal PSD. The inverse of the frequency at maximum residual PSD was used to 
obtain the dominant period. Second, the period at maximum peak-to-peak amplitude was tabulated by 
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measuring the zero crossings of the stacked waveform at each station (cf. ReVelle, 1997) in the same 
bandpass. The periods obtained using these two techniques show a very strong 1:1 correlation (Figure 
6), indicating that this methodology is robust in itself. 
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Figure 6. The dominant period correlation using two methods: PSD (vertical axis) and zero-crossings (horizontal 
axis). 


4 Estimating the Source Energy 

There are several empirical relations, relying on either the period at maximum amplitude or range and 
signal amplitude, which can be utilized in estimating source energy from infrasound measurements 
(Edwards et al., 2006). The yield estimates based on infrasonic amplitude are very uncertain in this 
instance as the propagation distances are much larger than is typical and outside the range limits where 
such relations have been developed (Edwards et al, 2006). In general, infrasonic period is less modified 
during propagation than amplitude (cf. Mutschlecner et al., 1999; ReVelle 1997; ReVelle 1974) and thus 
the period relationship is expected to be more robust. The Air Force Technical Application Centre 
(AFTAC) period-yield relations which are commonly used for large atmospheric explosions, are given 
by ReVelle [1997], as: 
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log(£/2) = 3.341og(P)-2.58 


E!2<\mkt 


( 1 ) 


log(F72) = 4.141og(P) - 3.61 E/2>40kt (2) 

Here, E is the total energy of the event (in kilotons of TNT), P is the period (in seconds) at maximum 
amplitude of the wavefonn. Since these relations were originally derived from nuclear explosions, the 
factor Vi must be incorporated in order to account for energy loss due to radiation for low altitude 
nuclear airbursts (Glasstone and Dolan, 1977). Even though there are a number of effects that may 
adversely influence and change the period at maximum amplitude during long range propagation of 
infrasound, this approach remains more robust than the maximum amplitude based relations, since it 
shows better agreement with energy estimates for bolide events which had their energies estimated by 
other methods (Silber et al, 2009; Brown et ah, 2002). 


5 Results and Discussion 

There are total of 17 detections, 16 obtained with PMCC and 15 obtained with MatSeis (Table 1). These 
detections overlap, except for the signal detected via MatSeis at the Bolivian station (I08BO), 17 509 km 
from the source. This signal, though very weak and short in duration (~30 seconds) compared to other 
signals (>185 seconds), shows a strong correlation to the bolide. The correlation indicators are the 
arrival time, the signal velocity, the dominant period and the apparent agreement between the observed 
and expected azimuth. The first arrival was detected almost two hours after the event at the closest IMS 
station, I39PW, at 04:39:51 UTC, while it took nearly 15 hours for the last bits of the signal to arrive to 
I08BO. Duration of the signal at each station (not including I08BO) was quite significant, ranging from 
3 minutes up to 41 minutes. All infrasound signals from the event show similar characteristics, such as 
long period and very low frequency content, consistent with a large blast radius and consequently a large 
energy source (ReVelle, 1976). Furthermore, average signal celerities are between 270 m/s and 320 m/s, 
indicative of stratospheric duct signal returns. 

The presence of high altitude winds affects the propagation of the signal in such way that it 
amplifies the downwind propagation, while it attenuates upwind propagation (c.f. Mutschlecner and 
Whitaker, 2010; Davidson and Whitaker, 1992; Reed 1969a). Most of the detecting stations are located 
east from the source and in October the stratospheric winds are predominantly westerly in the northern 
hemisphere (Webb, 1966). Average signal celerities (defined by the ratio between the horizontal 
propagation range and the travel time) are between 0.27 and 0.32 km/s, which is consistent with 
stratospheric duct signal returns. We also searched for possible antipodal signals, but found none. 

The geolocation ellipse (Figure 5), computed using azimuths and arrival times, points to 4.9°S 
and 122. 0°E with mean residuals of 2.9°. The source time estimated from this location is 02:52:22 with 
a residual of 1320 s. The accuracy of the source location strongly depends on the atmospheric wind and 
temperature profiles at the place and time of the event. 

To establish the best possible energy estimate of the Indonesian bolide, the average global period 
as well as individual periods, using both previously described zero-crossings and PSD methods, for each 
station were utilized. Table 2 shows the summary of energy estimates. The combined average periods of 
all phase-aligned stacked waveforms at each station produce a global average of 14.8 seconds (zero 
crossings method ) and 15.3 seconds (PSD method), corresponding to a mean source energy of 42.7 kt 
of TNT and 47.3 kt of TNT, respectively. Using the measurements from nine stations with the highest 
signal-to-noise ratio energy yield is 66.1 kt of TNT (zero crossings method) and 78.1 kt of TNT (PSD 


263 


method). The standard deviation of energy measurements across all stations is approaching the 
measurement itself, but this is expected because the signal usually emanates from different portions of 
the bolide trail as observed at different stations. Our best source energy estimate is 70 ± 20 kt TNT, with 
the error bounds representing the spread in the average from the different approaches (Table 2). 


Table 2. List of all detecting stations and their periods measured via two methods (zero-crossings at maximum 
amplitude in time domain and frequency at maximum PSD in frequency domain), as well as energy measurements 
for each station, where appropriate AFTAC relations were used (equation (1) or equation (2)). 


Energy estimate as a function of period I Energy estimate as a function of SNR 
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7.11 
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4.05 
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IS 18 

21.81 
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25.60 
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72.10 
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Energy estimate based on averaged global period 






Total 

14.81 

42.73 

15.27 

47.30 







6 Conclusions 

The Indonesian bolide of 8 October, 2009, detected infrasonically on a global scale, was perhaps the 
most energetic event since the bolide of 1 February, 1994 (McCord etal., 1995) and may have ex- 
ceeded it in total energy. We have no other instrumental records of this event other than casual video 
records of the dust trail emphasizing again the value of infrasonic monitoring of atmospheric explosive 
sources. Low frequency waves were observed at 17 IMS stations of the CTBTO network, making it one 
of the best infrasonically documented events (DTRA Verification Database, available at: 
http://www.rdss.info). 

Using an average impact velocity for Near Earth Objects (NEO) of 20.3 km/s, the energy limits 
(50-90 kt of TNT) suggested by this analysis correspond to an object 8-10 m in diameter. Given our 
upper limit in energy and a lowest possible entry velocity of 1 1 .2 km/s, the upper limit to the mass for 
this meteoroid is < 6000 tonnes. Based on the flux rate from Silber et al. (2009), such objects are 
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expected to impact the Earth on average every 10-22 years. Additional instrumental records of this 
unique event would prove valuable in understanding in more detail its interaction with the atmosphere 
and documenting possible local atmospheric perturbations. 

Additional instrumental records of this exceptional event, such as seismic, ground video 
recordings, satellite and possible meteorites, would prove valuable in understanding such occurrences 
and documenting possible local atmospheric perturbations. Since events like this one are rather rare, it is 
essential to maximize all aspects of such observations in order to validate propagation models at global 
scale, implement and better understand the spatial and temporal influences of atmospheric dynamics 
over propagation times, especially over long distances, and to evaluate energy yield formula and 
establish what information, not available via other techniques, can be derived from infrasonic 
measurements. 
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